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I. INTRODUCTION 

111 the past few years there has been increasing interest in the holographic duahty relating fluid dynamics and 
gravity. The link was first introduced in the 1970's with the development of the membrane paradigm approach [1], 
where the dissipative black hole horizon dynamics is recognized to closely resemble that of a viscous fluid. The 
connection has been made more concrete in the context of the AdS/CFT correspondence, where quantum gravity in 
an asymptotically anti-de Sitter (AdS) spacetime is shown to be dual to a certain gauge theory in flat spacetime in 
one lower dimension [2] . The gauge theory can be thought of as living on the timelike AdS boundary, in which the 
bulk spacetime is holographically encoded. As a consequence, the relativistic hydrodynamics of the gauge theory can 
be effectively described by the long time, long wavelength dynamics of a black hole in AdS [3] . The relativistic Navier- 
Stokes equations turn out to be equivalent to the subset of the General Relativity (GR) field equations called the 
momentum constraints, which constrain "initial" data on the timelike AdS boundary. Moreover, the incompressible 
Navier-Stokes equations describing ordinary, everyday fluids can be obtained by taking a particular non-relativistic 
limit of these results [4] . 

A key step toward a deeper understanding of the fluid/gravity correspondence can be found in the question of 
whether an asymptotically AdS spacetime is actually a required ingredient. Indeed, there have been a number of 
hints indicating that this is not the case. For example, the momentum constraint equations are not affected by the 
value of the cosmological constant, which suggests that the full asymptotic structure of the spacetime is un-important. 
Secondly, in cases where one can perform a hydrodynamic (long time, long wavelength) expansion of the equations 
describing the horizon dynamics, one finds they also have the form of Navier-Stokes equations [5-7]. Interestingly, 
one example where such an expansion exists is for a Rindler acceleration horizon in fiat spacetime [5, 8]. 

In a recent paper [9], a novel formalism was introduced to describe a holographic fluid theory dcflncd on an 
arbitrary timelikc surface in a general spacetime with a causal horizon. On this surface, one flxcs the boundary 
condition that the induced metric is flat, and in the spirit of the Wilsoiiian approach to the renormalization, the 
asymptotic physics outside this surface plays no role. Moving this surface between the horizon and the asymptotic 
boundary can be thought of as a renormalization group flow between a boundary fiuid and a horizon fluid. In [10] 
the authors considered the specific case of perturbations about a Rindler metric, taking the timelike surface to be one 
of the family of hyperbolas associated with the worldlincs of an accelerated observer. Working in the non-relativistic 
hydrodynamic expansion, the authors presented a geometry that is a solution to the Einstein equations if the data 
on surfaces of Tc satisfy the incompressible Navier-Stokes equations. Alternatively, one can consider the physically 
inequivalent near-horizon expansion in small Tc and obtain the same results. 

Beyond the connection between the classical Navier-Stokes equations and a classical geometry, these works actually 
suggested the possibility of an underlying holographic duality relating a theory on fixed Vc to the interior bulk of the 
Rindler spacetime. A first step toward a detailed study of the behavior of this dual system was taken in [11], with 
the introduction of an algorithm for constructing the geometry and the explicit expression for the viscous transport 
coefficients to second order in the hydrodynamic expansion. 

In this paper, our main goal is to probe further the dual theory by asking what effect higher curvature terms in the 
dual gravitational theory have on the transport coefficients of the fiuid dual to the Rindler geometry. In the AdS/CFT 
correspondence, such terms are associated with quantum corrections or other deformations, which modify the values 
of the transport coefficients. Remarkably, we show here that the shear viscosity of the Rindler fiuid is not modified 
if higher curvature terms are introduced. Equivalently, at lowest orders in the non-relativistic expansion, the dual 
metric solution has the property of being a solution to GR and to any higher curvature theory of gravity. The first 
place the higher curvature corrections appear is in the second order transport coefficients of the fiuid. Working in the 
case where the higher curvature theory is Einstcin-Gauss-Bonnct gravity, wc calculate some of these coefficients. 

The plan of this paper is as follows. In Section II, we describe the general construction of the solutions developed 
in [9-11]. In Section III, wc explicitly show that the shear viscosity of the dual theory is unchanged when generic 
higher curvature terms arc added into the gravitational action and wc discuss the differences between this calculation 
and previous literature on the AdS/CFT correspondence. Section IV is devoted to the calculation of the second 
order transport coefficients in Einstein-Gauss-Bonnet. We then conclude with a discussion of the implications of these 
results and their possible connection to approaches using the local Rindler geometry as a tool for a thermodyiiamical 
derivation of gravitational dynamics. 

II. GENERAL SETUP 

Wc want to construct a Lorentzian geometry that acts as the holographic dual description of a fiuid fiow in d + 1 
dimensions. Based on the holographic principle, we expect the fluid is defined on a d + 1 dimensional timelikc surface 
Sc embedded in a d -I- 2 dimensional bulk spacetime. We choose the timelike surface to be defined by fixed bulk radial 
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coordinate, r = Tc- Wc also specialize to the case where the fluid moves on a flat background. In this case, the induced 
metric on Sc should be flat as well, e.g. 

j^^dx^'dx" = -$(r,)dt2 + e2*('^<^)da;,dx*, (1) 

where $ and 5* are some functions of r. We use the notation that coordinates on the hypersurface Sc are = {t, cc'), 
where i = l...d. The {d + 2) dimensional bulk coordinates are defined with the notation x^ = (t,x%r). The final 
requirement is that the bulk spacetime must contain a regular, stationary causal horizon. The bulk spacetime therefore 
has a timelike Killing vector field, which becomes null on the horizon. The full bulk metric therefore has the general 
form [9], 

ds"^ = -$(r)dt2 + 2dtdr + e2*('-)dx»da:*, . (2) 

where at some radius r = rh there is a horizon where $(r) ~ and the timelike Killing vector = (f^t)^ becomes 
null. If one considers quantum field theory on the background (2), one finds equilibrium thermal states associated 
with the presence of the horizon. For example, one can compute the Hawking temperature (in units where h ^ c ~ 1) 



where the surface gravity k can be defined via bX^ ~ i^X^ ■ Dividing by the redshift factor at rc, \J—gtt = \f^Wc) 
yields the local Tolman temperature 

^loc ~ - — 7==- [V 

There is also an associated Bekenstein-Hawking entropy proportional to the cross-sectional area of the horizon 

5BH = 4^e'^*('-'>), (5) 

where here and throughout we use units such that IGttG = 1. We want to identify these thermodynamical properties 
with the thermodynamical properties of the dual fluid in d + 1 dimensions. Therefore, the general metric can be 
thought of as the dual geometrical description of an equilibrium thermal state associated with some lower dimensional 
theory defined on the surface r = Vc- 

The metric (2) can describe many different black hole solutions. Here we will focus on the special case of a region 
of flat {d + 2) dimensional Minkowski spacetime in "ingoing Rindler" coordinates 

ds^ = -rdt'^ + 2dtdr + dxtdx\ (6) 

where in terms of the above parametrization, $(r) = r and \E'(r) = 0. The null surface r = acts as a horizon to 
accelerated observers, whose worldlines correspond to surfaces of constant r = re- 

Although the Rindler metric is just a patch of flat spacetime, the associated quantum field theory on this background 
has many of the same properties as a black hole solution. In particular, surfaces of r = rc have a local Unruh 
temperature 



Strictly speaking, a Rindler horizon does not have a Bekenstein-Hawking entropy density. However, one can assign 
the Rindler horizon this entropy based on the holographic principle, or, more concretely, take the entropy to be the 
thermal entanglement entropy of the quantum fields in Rindler wedge [15]. This statistical entropy scales like an area, 
but is a UV divergent quantity. If a Planck scale cutoff is chosen appropriately, the entanglement entropy agrees with 
the Bekenstein-Hawking formula, i.e. 

s = 47r. (8) 

Given the existence of an equilibrium Unruh temperature and a Bekenstein-Hawking entropy density, the metric (6) 
can be thought of as a dual geometrical description of a perfect fiuid in one lower dimension. This duality can be 
formalized by considering the Brown- York stress energy tensor [12], which in GR takes the form, 

T/i^J" = 2{Kj^,~K^,), (9) 

where Kf^i, = ■^CnI^v and Cm is the Lie derivative along the normal to the slice N^. One can show that (and 
its generalization for higher curvature gravity) is indeed equivalent to the stress energy tensor of the perfect fiuid with 
a rest frame energy density p and pressure P. In this case 

P = 0, P=^. (10) 
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III. EQUIVALENCE OF VISCOUS HYDRODYNAMICS IN EINSTEIN AND HIGHER CURVATURE 

GRAVITIES 

A. The seed metric 

In this section we will argue that the first order viscous hydrodynamics of the fluid defined on Sc is independent 
of whether the dual gravitational theory is Einstein or some higher curvature generalization. In order to study the 
hydrodynamics of this fluid, we must perturb the background Rindler geometry. To start, we review the formalism 
for perturbing the Rindler metric developed in [11]. The first step is to make a set of coordinate transformations 
to obtain a new metric (or class of metrics). These transformations should keep the induced metric at Tc flat. The 
transformed metric should also preserve a perfect fluid form of the stress energy tensor associated to the slice, as 
well as the time-like Killing vector and the homogeneity in the direction. It was shown in [11] that these set of 
conditions uniquely identify the two diffcomorphisms, namely a boost and the translation. 

The boost of the metric takes the form, 

y^ct^V^ct-7l3,x\ ^x' -^P'^ct+{l-l)^x\ (11) 

where 7 = (1 — /?^)^^/^ and Pi = r^^^^Wj is the boost parameter. The linear shift of the radial coordinate and 
re-scaling of t, which moves the horizon from r = to an r = r/i < Tc, is instead 

r^r-Th, {l-rh/r,)-^'h. (12) 

The resulting metric for the flat spacetime is 

ds^ ~ ^ — \ J — I -\ ^ — dtdr ^ ^ — d/T^dr 



dx'dt + S^j - ' , . '— dx'dx^. (13) 



We now want to investigate the hydrodynamic system dual to the above metric. To do that, we need to consider the 
dynamics of the metric perturbations within a hydrodynamic limit. One can perturb (13) by promoting the spatial 
velocity and horizon radius to be functions of space and time: v'^{t, x*) and rh{t, x'). Now the metric is no longer flat 
and no longer a solution of the vacuum Einstein equation. However, one can introduce a particular non-relativistic 
hydrodynamical expansion [4, 13] in terms of a small parameter e, 

v' ev'{tx\eH) P e'^P{ex\eH), (14) 

where the non-relativistic pressure P(t,x^) is defined in the following way as a small perturbation of the horizon 
radius, ^ 

rh=0 + 2P + O{e^). (15) 

Using (14) one scales down the amplitudes (e can be thought of as the inverse of the speed of light), while at the 
same time scaling to large times t and spatial distances x^ . This corresponds to looking at small perturbations in the 
hydrodynamic limit. 

Expanding the metric (13) out to O(e^) in this manner yields the "seed metric" solution originally found by 
Bredberg, Keeler, Lysov and Strominger in [10], 

ds^ = -rdt^ + 2dtdr + dxidx"- 

2(1-—^ v,dx'dt - —dx'dr 

1 - ^ ) [(^2 + 2P)dt^ + ^dx'dx^ +(— + --) dtdr. (16) 



/«2 


2P\ 







^ Note that the e expansion is performed in such a way that at zeroth order = rj, = so that the standard Rindler metric (6) is 
recovered. Also, there is no scaling of bulk radial derivatives. 
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The seed metric is the unique singularity-free solution to the vacuum Einstein equations up to 0{e^), provided div'^ = 0. 
As required, the induced metric on the slice r = is flat. 

In GR, the momentum constraint equations on the surface Sc can be expressed in terms of the Brown- York stress 
tensor 



R.aN'' = d^T^J = 0. (17) 
At second and third order in e, momentum constraint equations are 

R^fN^ = r- V2^(2,3) ^ ^1/242.3) ^ 0, (18) 

while the Brown- York stress-tensor for the seed metric is given by [10] 



+ v,vj - 2rcdiVj dx'dx^ + 0{e^) . (19) 

Then, at second order, using the expression in (19), the momentum constraint equations (17) reduce to the incom- 
pressibility condition div' = we discussed above. At third order one finds the Navier-Stokes equations with a 
particular kinematic viscosity 

dm + v'djv, + d^P - T.d'^v, = 0. (20) 

Therefore, imposing the the incompressible Navier-Stokes equations on the fiuid variables guarantees the dual metric 
is a solution to the field equations. 

Noticeably, these results can be obtained as a non-relativistic expansion of a relativistic viscous fluid stress tensor. 
To see this, we work in the relativistic hydrodynamic expansion in derivatives of the fluid velocity and pressure: du 
and dp. Then, at first order, the relativistic viscous fluid stress tensor has the form. 

Til"" = puy^u, + p/v. - 277/v, - ih^,,{^xu^)■ (21) 

Here /i^i/ = 7^,^ -I- u^Ui,, while K^j^^ = h^h1d(^\u„) is the fluid shear, r] the shear viscosity, and ^ the bulk viscosity. 

The viscous terms above are written in the Landau or transverse frame [17], which can be defined as a condition 
on the first order part of the stress tensor 

yfluid (1)^^ (22) 

This frame is constructed so that the viscous fluid velocity is deflned as the velocity of energy transport. The seed 
stress tensor in (19) follows from the e expansion of (21), if we identify 

^ -Xrcy), p = + O(e3), p=^ + JL, ^ = 1. (23) 



This is consistent with the earlier equilibrium calculation of p and p in (10). Note also that the bulk viscosity term 
in (21) actually drops out and bulk viscosity is not an independent transport coefficient. This is due to the fact that 
at viscous order we can impose the ideal order equation ~ 0, which follows from p = and continuity. 



B. Higher Curvature Gravity 

Now we want to study how the hydrodynamics of the fluid is modified when the gravity theory is not GR, but 
instead some theory with higher curvature terms. The first question is whether we need a new. modified seed metric 
in a higher curvature theory of gravity. Interestingly, we can show that the seed metric (16) and its 0{t^) correction 
is a solution to a wide class of higher curvature gravity theories at lowest orders in the e expansion. 

We start by noting that the flat, equilibrium Rindler metric at zeroth order is a vacuum solution to both Einstein 
and higher curvature gravity theories. The higher curvature terms could be thought of as modified gravity theories 
in their own right or they can be seen as quantum corrections to Einstein gravity in an efiiective field theory picture. 
Here we will not consider exotic theories involving inverse powers of curvature invariants. 

As a first example of a higher curvature theory we consider Einstein-Gauss-Bonnet gravity (in the absence of a 
cosmological constant), defined by the action 

las = / d'^+^x^g [i? + a (i?2 - 4i?cDi?^^ + i^CDBFi?^^^^)] , (24) 
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where a is the Gauss-Bonnet couphng constant. We consider d > 3 since for d < 3 the Gauss-Bonnet term is 
topological and does not affect the field equations. The interest in looking at a Gauss-Bonnet term is twofold. Such 
a term arises in the low energy limit of string theories. Secondly, Einstein-Gauss-Bonnet gravity is notable because 
even though the action is higher order in the curvature, for the unique combination of curvature invariants in the 
second term of (24), the field equations remain second order in derivatives of the metric. 
Varying this action with respect to the metric yields the field equations. 

Gab + 2aHAB = 0, (25) 

where the Lovelock tensor Hab is 

Hab = RRab - 2RacRb - 2R^"Racbd + Ra^^'^'Rbcde - ^gAB (R^ - ^RcdR^'' + RcdefR^'''"') ■ (26) 

Now, using the seed metric, the first non-zero components of the Riemann tensor Rabc^ a-rc at O(e^). If we examine 
the Lovelock tensor, (26), it is clear that the first contributions from the Gauss-Bonnet terms can only appear at 
O(e^) at the lowest. A similar conclusion obviously holds for Lovelock gravities [14], which are the extension of the 
action (24) including contributions with higher powers of the curvature but still yielding 2nd order field equations. 

The field equations of other higher curvature theories of gravity generally involve covariant derivatives of the 
Riemann tensor and its contractions. These are no longer second order in metric derivatives. At second order in the 
curvature the gravitational action has the form 

7AB : on nABCD\ 



PiR' + P2RABR'''' + P^RABCoR^'^'^n ■ (27) 
The field equations can be expressed in the form Gab = ^abj where 

Sab = Pi (^RRab - ^a^bR + 9AB{nR - ^R'] 

P2 {9abRcdR^'' + ^S/c^bR'a - 2nR^, - qabOR - AR^^Rcb) 

+ {qabRabcdR^''^'' - ^RacdeRb^'''' - SDRab + aR + SR'IRcb - SR^'^Racbd) ■ (28) 

Let's consider the possible terms that can appear at the lowest orders in e. First, the second covariant derivative 
terms of R could in principle contribute Pi corrections at 0{e^). However, the Ricci scalar R = g^^ Rab can be 
expanded out as follows, 

R = + 2g^'Rtr + 2g''Ru + 2g"R„ + .g^'^i?„ + . (29) 

Before imposing incompressibility, one can show that the only non-zero component of Rab at O(e^) is 

Rtt = \d,v\ (30) 

However, for the background Rindler metric (6), the zeroth order <?**(o) is zero, so the Ricci scalar R is in fact higher 
order. Since one cannot form a scalar constructed from P, dt, and di with odd powers of e, we expect R is of 
O(e^). For instance, the spatial vector Ru is O(e^), but this multiplies 5**, which is 0(e). Therefore, R is 0{t^) and 
its covariant derivatives are of the same order or higher. 

The remaining terms of interest are the 'CIRab and VcVsi?^ terms proportional to P2 and P-^. We know that Rab 
a priori has non-zero components at O(e^) and 0{t^). The question is whether the radial derivatives and background 
connection for the Rindler metric (6) allow the above two terms to also contribute at these orders in e thereby affecting 
the hydrodynamics at these orders. This we checked with an explicit calculation. The result is again negative. 

Thus, as a general principle, higher curvature corrections to the Einstein equations come in at O(e^), at least 
when we perturb the fluid dual to the fiat Rindler spacetime geometry. Terms of even higher order in the curvature 
(schematically ^ i?", where n > 2) will typically appear at even higher orders. This includes the often studied case 
of f{R) theories, when / can be expanded around the Hilbert term: f = R + R^ + R'^ + ■ ■ ■ . 

As a result, the solution to the higher curvature theories at the lowest orders O(e^) and O(e^) is the same as 
the GR solution found previously [10, 11]. Since all the higher curvature quantities vanish at the lowest orders, 
this solution has the property of being approximately strongly universal [16] . The explicit solution at 0{e^) can 
be constructed from the algorithm for Einstein gravity given in [11], which we will expand upon and generalize to 
Einstein-Gauss-Bonnct in the next section. At the present, we note that the equivalence of the solutions to 0{e'^) 
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implies that the 1st order viscous hydrodynamics of the dual fluid is the same both in Einstein gravity and its higher 
curvature generalizations. In particular, the incompressible Navier-Stokes equations (20) are the same in any theory, 
with the kinematic viscosity fixed to be Tc. Furthermore, as Compere, et. al. pointed out, it is clear that the 
non-relativistic e expansion is capturing the non-relativistic limit of a relativistic fluid theory whose full structure is 
unknown. Nevertheless, the e expansion seems to be able to capture some of the transport properties of this fluid 
theory. In particular, the shear viscosity of the relativistic fluid, 77, is apparently flxed to be 1 (or (IGttG)^^ if we 
restore the gravitational constant). 

One may worry about using the non-relativistic limit to draw conclusions about the properties of the relativistic 
parent fluid. However, we can show that our analysis of higher curvature terms can be extended to the relativistic 
hydrodynamics. The flrst step is write the metric (13) in a manifestly boost covariant form. This metric turns out to 
be 



ds^ = —(1 + p^{r — rc))ui_iUudx^dx'^ — ipu^dx^dr + hf^^dx^dx'^ . (31) 
In this line element we have replaced with the relativistic pressure p using the general formula 

p^—^=. (32) 

V1"c - Th 

Expanding and p in terms of and P using (23), the metric (31) reproduces the seed metric up to O(e^). In 
addition, if we compute the Brown- York stress tensor at r = for this metric (31), we find directly 



T^ydx^dx^ = ph^^dx^dx'^ , (33) 

which is the ideal part of (21) with p = 0. 

To perturb in this case, we now treat u'^{x^) and p{x'^), but leave fixed. The metric is no longer a solution to 
the vacuum Einstein equations, but one can expand and work order by order in derivatives of u'^ and p as discussed 
earlier. This follows the standard approach used in the fluid-gravity correspondence [3]. 

We now have that (31) is a zeroth order solution, i.e. Rab =0-1- 0(A), where the parameter A counts derivatives 
of u'^ and p. Therefore, Rabc^ ^ and the curvature squared terms in (28) must appear at O(A^). The other 

terms involve the covariant derivatives of the Ricci scalar and tensor. The generalization of (29) is 

= -f 2ff'^^i?W + gM^i^W . (34) 

From (31) we find 




d[f_,pu^) + Dp u^Ui, + ]^p{d\u^)u^u^ + pui^f^a„) , (35) 

where we have defined D = m^'S^ and ^ u^dxu^^. Since g^"^ = /i^"^, which projects orthogonal to u^^, R^^^ = 0. 
Finally, the fact that the remaining terms URab and '^c^bRa are also of O(A^) can be shown by explicit calculation 
as before. 

Therefore, we conclude again that the higher curvature terms affect only the second order viscous hydrodynamics. 
The equilbrium stress tensor will be given by Eqn. (33) in any higher curvature theory of gravity. This follows just 
from the fact that the zeroth order metric (31) is a solution in any theory. Computing the 0(A) corrections to (31) 
and (33) confirms that 77 = 1 and the bulk viscosity is not a transport coefficient, but we will save the details for 
another paper [18]. 

In higher curvature theories, the entropy is given by the Wald formula [19]. In general, Bekenstcin-Hawking area 
entropy will be modified by the higher curvature terms, leading to an expression that can depend on both the intrinsic 
and extrinsic geometry of horizon. However, since we are working with a Rindler horizon in flat spacctimc, all these 
corrections vanish and the equilibrium entropy density s remains 47r. The ratio ry/s = l/47r was first derived in the 
context of the AdS/CFT correspondence [20]. It was shown that the ratio goes to this value for any infinitely strongly 
coupled holographic gauge theory fluid with an Einstein gravity dual [21]. On the gauge theory side, the number of 
colors iV — 00 and the 't Hooft coupling A — > 00. This is essentially a classical limit; quantum corrections to the rj/s 
ratio at finite N and A, which can be calculated in specific string theory realizations [22], correspond to specific higher 
derivative corrections to the dual gravitational theory. Another approach is to work outside the context of particular 
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string theories and consider a generic higher curvature gravity action of the form given in (27). In this case, it has 
been shown [23, 24] that ratio changes to 

1 (l_8/?3). (36) 
s 47r 

This result holds in five spacetime dimensions and to linear order in the Pi, which are effectively suppressed by powers 
of the Planck length. It is also important to note that while the ratio is unchanged when (3^ = 0, both rj and s do 
depend on /3i^2- Finally, in the special case of Einstein-Gauss-Bonnet, (24), /33 = a. Given the nice properties of 
this theory (linked to the field equations remaining 2nd order in derivatives), one can work non-pcrturbatively and 
consider finite a corrections which allow the ratio to approach zero. 

It is then remarkable that in the case of a flat Rindler background there is no higher curvature correction to the 
ratio or to the viscosity itself. The viscosity is protected against quantum corrections or other deformations to the 
dual theory. At a technical level, the difference is that the result (36) follows by considering perturbations around a 
background asymptotically AdS black brane solution in the higher curvature gravity theory. In Einstein-Gauss-Bonnet 
gravity with negative cosmological constant, this solution is [25] 

ds^ ^ N^f(r) - -i-rfr^ + r^dx^dx' (37) 
/('') 

where N is some constant and 



^(") = 4^l^^"V^"^''v ^^^^ 

with r/i the value of the horizon radius. In this solution, thermodynamic quantities such as the temperature and 
entropy density depend explicitly on a, which translates into the calculations of the entropy and shear viscosity. 

In contrast, in the Rindler case the metric does not depend on a and the Unruh temperature and entanglement 
entropy are kinematical quantities in the sense that they are independent of the underlying gravitational theory. The 
shear viscosity seems to have the same behavior since it is also unaffected by the choice of gravitational dynamics. 
This is further evidence for the picture of ij/s = 1/Att as a kinematical property associated with entanglement in 
Rindler spacetime [26]. 

IV. SECOND ORDER TRANSPORT COEFFICIENTS 

Now let's consider the hydrodynamic expansion at higher order in derivatives. Here we expect the gravitational 
dynamics to affect the hydrodynamics of the dual fluid. To second order, O(A^), the general stress tensor for a 
relativistic fluid with zero energy density (hence incompressible) has the form [11] 

+ csa^^ D Inp + c^Dj^ In p Inp, (39) 

where D ~ -u^S^, Dj^ = P^d^, and VL^u — P^P^d^xu^y There are also viscous corrections to the energy density p at 
this order, which can be parameterized as 

P = hiK^^K'"' + b2n^M^" + hDlnpDlnp + b^D'^ \np + b^Dj^ InpO^i" \np. (40) 

The Ci, i = 1..6, and bj, j — 1..5, are the possible new transport coefficients. When one expands these expressions in 
powers of e, many of the second order transport coefficients appear at O(e^) in a general non-relativistic fluid stress 
tensor. 



T^^f'^i ^^Ux^dx" = 



b b 



df 



(v^ + P) + 2?7ref(,a,)P + c^rl/^d^djP + ^r^/V.^a^- 



, '^3 3/2 k "^2 3/2 r, 2 c,2 a 2 



I c 2 



dx'dx^. (41) 
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Here aij = 2d(^iVj) and = 2d[iVj]. Only C5, Cq and 63, 64, and 65 are absent at this order in the e expansion. 

We argued that O(e^) is the first to receive corrections from any higher curvature terms in the gravity theory. In 
the next section, we wih solve for the fourth order (non-relativistic) metric in five dimensional Einstein-Gauss-Bonnet 
gravity. With this result in hand, wc will use the corresponding Brown- York stress tensor to read-off various second 
order transport coefficients for the dual fluid. 



A. Constructing the Einstein-Gauss-Bonnet solution 



We first outline the construction due to [11], where one starts with the metric solution at 0(e" In practice, 
the first n is 3 , i.e. one starts the process with the seed metric solution (16). Wc then want to add to the metric a 

new piece g^^g that solves the field equations to 0(e"+^). Since radial derivatives carry no powers of e, the addition 

(n) 

of produces a change in the bulk curvature tensors at the same order. This is effectively a perturbation around 
the zeroth order background Rindler metric (6). We work in the gauge where 

t2=0, (42) 

for all the contributions with n > 3. With this choice, we find that changes in the Einstein tensor 6Gab = SRab — 
■^g^B^R have the form 



5G^^ = -r5G':i> = {2rdlg^ + 9.5!"^) • (43) 

We define (?-"^ = 6^^ g^^^ and (5G-"^ = d^^SG["\ In contrast, there is no change to the Lovelock tensor (26) at the 
same order n since the curvature of the Rindler background is zero and any term in the variation would contain some 
factor of curvature at zero order. 

We want to find the that cancels out the 0(e") part of the field equations arising from the pre-existing solution 
at 0(e"^^). That is, we require 

SG^^l + G^;i+2aH2i=0 (44) 

where the hat denotes the parts of the curvature arising from the pre-existing solution. In order for this set of 
equations to be consistent, one must impose the integrability conditions 

mUrEl:^=0 (45) 
E^^+rElf=0 (46) 
driMr'^ + rEl'^'^ ) + (l/2)-B(;') = (47) 

where we have defined e'x% ^ G% + 2aH'xB- These are consistent with the Bianchi identity and (45) follows from 
the conservation of the Brown- York stress tensor extended to Gauss-Bonnet gravity [27], i.e. 

{Gau + 2aHA.)N^ = d^T,,„ = 0, (48) 

where 

T^, = 2(i^7^, - /V) + 4a(J7^, - 3J^, - 2P^p,,KP^). (49) 

The symbol P^pua — Rp.pua + 2Rp\^i,^„^^ — 2Rp'^^^„^p + R'y^[iy^a]p is the divergence free part of the induced Riemann 
tensor and can be neglected here because we work with a fiat induced metric, while 

V = I {2KK^,K: + K^xK^^K^, - 2K^,K"^Kxu - K^K^,). (50) 
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Using (43), one can solve the differential equations subject to two conditions: (i) that g^"^ = at r = Tc (the metric 
on Sc remains flat) and (ii) that there is no singularity at r = 0. The resulting solution is 



(") 
9tt 



(1 - r/r,)^^i") {t, x^) + j'^ dr' j ' ^' dr''^ (^"^ " 44<r^ - 2r^(';)) (51) 
ffi"^ = (1 - rlr^)Ft~^ it, x^)-2 dr' dr"^"' (52) 



9^' ^-2 dr''- dr"(E!(;^+2a{H^-^Ht;})] (53) 



where Ff"\t,x'') and F^"'\t,x''') are arbitrary functions. 

These two remaining functions can be fixed by imposing gauge choices on the Brown- York stress tensor of the fluid 
(49). The addition of the new metric piece at 0(e") has the following effect on the extrinsic curvature at the same 
order 



(54) 



implying that 



By explicit calculation, we verified that there is no corresponding 0(e") variation of the J^^ part of the stress tensor. 
Thus, the variation STjfl^ comes only from the linear part in the extrinsic curvature: 



= S.: ( ^^^^ + ^ /" dr'iRt"^ + ^Hif)] - ^ dr'(i?(;) + 2a7?(")). (56) 



rc V ' c Jo 

The complete stress-tensor has the form 

r^i:^ = ST^:^ + 2 (k^^h^u - Kil^) + 4a (i(")7^, - 3j(^)) , (57) 

where as before, the hat notation indicates the part of the stress-tensor originating from the solution at 0(e"^^). 
The function F^"'\t,x^) is fixed by imposing the Landau gauge condition (22) order by order in the non-relativistic 
expansion. This plays a role only at odd orders in e. The other function Ff"\t, x"^), which appears at even orders, is 
fixed by requiring that there are no higher order corrections to the definition of the non-relativistic pressure, i.e. the 
isotropic part of is 

at all orders. 



B. Solution to 0{e^) 

We now apply the algorithm to solve for the metric to 0{e^). One first starts with the seed metric solution (16) and 
constructs the solution at O(e^). As we argued earlier, the corrections due to the Gauss-Bonnet coupling constant 
arise at 0{e^). Therefore, the Gauss-Bonnet terms do not contribute and the solution reduces to the GR one found 
previously in [11], where the only no n- vanishing component is 



(3) 

9u 



2r, 



2u 



2P) — + Ad^P - (r + rc)d^v^ 



(59) 
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The next step is to compute the R^^g and H^^g using this metric. Via direct calculation of the Lovelock tensor 



(26), we find that 



(60) 



with all other components of H'^g equal to zero. At even order in e. Ru ~ and as a result g^^^ = 0. The remaining 
components to compute are R^^^ , rI^ , , and i?^-^'' , which we will not display explicitly here. 
Using (60), the solution for gl-f^ in (51) reduces to 



(4)/ 



and we find that 



g[t^ = {l-r/n)Fr{t,x^) + 



(4)/ 



8r, 



(r — rc)^(r — + 2a) 



-UJklUJ 



The gauge condition on the stress tensor (58) fixes 



q c: p2 

^ -.4 , ^ py2^^_ 2r^y^d^v' 



8r 



2r, 



(61) 



(62) 



(63) 



Note that in these expressions we have imposed incompressibility diV^ = 0, used the Navier-Stokes equation (20) to 
eliminate time derivatives of Vi, and imposed 



which also follows from (the divergence of) Navier-Stokes. Meanwhile Eqn. (53), yields 



(64) 



^v,vj{v' + 2P) + -v^^^,)P - Ad,djP - ^<J^k<j\ 



r — 5rc + 12a j, 
Wifco; j 



+ '^fe(''^j) ^^'^ + — 4 — 

1 2 1 / 2 \ "-"^ - 

- —v(,d^^v ~7rr<yij{v +2P) + —d.,jujkii- 



We now use (57) and (56) to find the stress tensor components T^P and T^^' . The non-zero components of the j/i 
tensor are 



(65) 

(4) 



1 



24r 



1/2 ' '^i] 



12r: 



3/2 



Using this result, we find 



(66) 



(4) _ -3/2 



v^{v^ + P)- -fa,ja^^ - r^aijV^v^ 



and 



i;,i;,(i;2 + P) + 2r,i;(,9,)P - 4r%djP ^^(^1 + ^^ a.fea^ - (^1 + 



3r r 

+ rla^diOj)^ - 2rlv(id'^Vj) + —^d'^o^ - rcVi^idj^v^ - -^o^v^ 



(67) 



(68) 



Note that the T^P has no a corrections. They cancel out and the energy density Tfj^^u'^u'^ is not affected by a at 
fourth order. Comparing with the general form of the fluid stress tensor (41) we read off that 



2a 



bi = -Vr'c, 62 = 0, ci = -2V^ 1 + _ 1 C3 = -4V^ 1 + _ 1 C2 ==C4 = -4^^ 



3a 



(69) 



as expected, there is no change in the value of 77 = 1. However, the Gauss-Bonnet term does modify the two transport 
coefficients ci and C3 from their purely GR values. 
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V. 



DISCUSSION 



We have argued that higher curvature corrections to the Einstein equations always come in at 0{e'^) in the non- 
relativistic hydrodynamic expansion and at O(A^) in the relativistic Knudsen number expansion, at least when we 
perturb the fluid dual to the flat Rindler spacetime geometry. Hence, the solution to the higher curvature theories 
at the lowest orders is the same as the GR solution found previously [10, 11]. Working in the specific case where 
the higher curvature theory is Einstein-Gauss-Bonnet gravity, we then showed explicitly that the 1st order viscous 
hydrodynamics of the dual fluid is the same both in Einstein gravity and its higher curvature generalization, while 
the effect of the higher curvature corrections shows up in the second order transport coefhcients of the fluid. We 
calculated some of these transport coefficients and found that two of them depend on the Gauss-Bonnet coupling 
constant. It would be interesting to complete the relativistic calculation outlined in Section III B in order to find all 
the second order transport coefficients in both the Einstein and Einstein-Gauss-Bonnet examples. 



The approximate strong universality [16] of the seed solution about which the hydrodynamic expansion is made is 
an interesting result. The lack of a higher curvature correction to the viscosity implies that it is protected against 
quantum corrections or other deformations to the dual theory. One way of thinking about these results is to note that 
shear viscosity and entropy density typically scale like T'^, where T is the equilibrium temperature of the thermal 
system. In AdS/CFT. this temperature is given by the Hawking temperature Th of the black branc solution, which 
would depend in this case on the Gauss-Bonnet coupling constant, due to the non-trivial curvature of the background 
solution. In contrast, the shear viscosity and entropy density are constants independent of the temperature in the 
Rindler case. This suggests the two types of holographic duality are different. 

The independence of the fluid/Rindler holographic duality from the asymptotic geometry makes this correspondence 
interesting beyond the AdS/CFT context. For example, the Rindler metric is associated to an accelerated observer 
in the locally flat surroundings of any point in spacetime. Therefore, one can ask whether the flat spacetime duality 
can be applied locally and then possibly used to patch together a holographic description of any spacetime [11]. 

In a similar manner, the local Rindler system is also crucial the idea that gravity may emerge from the holographic 
hydrodynamics of some microscopic, quantum system [8, 28-30]. Here one assumes that the local Minkowski vacuum 
state carries a finite area entanglement entropy which can be holographically identified with the entropy of the local 
Rindler horizon. Perturbations to the horizon system are assumed to obey an entropy balance law, relating a change 
in the entropy to the "heat" associated with a flux of matter, plus an internal entropy production term from shear 
viscosity. Demanding that this equation holds at each point in spacetime then yields the GR Einstein equation and 
fixes the shear viscosity to entropy density ratio to be l/in. Also, bulk viscosity does not appear as an independent 
transport coefficient [8, 31], which is strikingly similar to the viscous hydrodynamics of the global Rindler fluid. 

However, extensions of this type of derivation to /(i?) gravities [29-31] require the horizon entropy to depend on 
the curvature, which inevitably leads to the same behavior in the shear viscosity. Ultimately, while the metric around 
any point is flat, the curvature itself does not vanish at the point. This fact means that one cannot simply import 
results pertaining to perturbations of the globally flat Rindler solution into the locally flat patch. On the other hand, 
there is some evidence that the inherent fuzzincss in the local Killing vector, which is associated with the local notion 
of thermal equilbrium, may be of the same order of magnitude as higher curvature corrections [29]. If this is the 
case, the approximate notion of a local fluid would not be affected by these corrections. It would be interesting to 
investigate further the relationship between the fluid/Rindler correspondence and these ideas of emergent graviational 
dynamics. 

Finally, to conclude, we want to point out an interesting duality here between the relativistic A expansion in 
derivatives in the holographic theory and an effective field theory expansion of the bulk gravitational theory. First 
note that the A expansion is equivalent to an expansion in small dimensionless Knudsen number, which is defined as 
Kn — :^2£e.^ where €nifp is the mean free path associated with the microscopic system and L is the characteristic size 
of the perturbations to the system. Secondly, although the bulk gravity theory is non-renormalizable, it is still valid 
as an effective theory when the dimensionless ratio of the Planck length to the radius of curvature, , is small. 

The effective action is given as an expansion in this ratio. At zeroth order there is some cosmological constant, at 
first order, the Hilbert term, and then the pieces higher order in curvature invariants. 

In the duality, the scale of perturbations L in the system on r = Tc is linked to the scale -Rcurv of perturbations 
to the flat bulk spacetime. As we have seen, a flat spacetime is dual to the fluid in equilibrium, Einstein gravity 
dual to the viscous hydrodynamics characterized by a shear viscosity, and second order transport coefficients linked 
to curvature squared terms. It is tempting to associate the universality of the shear viscosity with the universality of 
the Hilbert action at low energies and take i'mfp ^ ^pianck, the scale at which gravity is strongly coupled. This line 
of reasoning also suggests it may be interesting to consider a seed metric constructed from the region of a de Sitter 
spacetime where there is also a causal "observer dependent" horizon and the associated thermodynamics. What effect 
does a non-zero cosmological constant have on the dual fluid? 
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